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ABSTRACT 



An analytical investigation was made of combined forced 
and free turbulent convection in a circular vertical tube with 
volume heat sources and constant wall heat addition. 

The IBM 1620 digital computers located at The George Wash- 
ington University and also at the United States Naval Academy 
were used to solve the basic equations in which the parameters 
are: Prandtl number Pr; Rayleigh number Ra; friction Reynolds 

number Re*; and the volume heat source parameter F. For fixed 
values of these parameters the solution gave the fully developed 
velocity profile, the temperature profile, and the pressure drop. 
The program also solved for the mixed-mean-to-wall temperature 
difference, Nusselt number, and Reynolds number. The following 
values of these parameters were investigated: Pr = 1, 10, 100; 

Ra = 0, 3^, 5**, 10^; Re* = 103, 10 2 *, 105; and F = .1, .5, 1, 10. 

The results of laminar flow problems were investigated to 
check the validity of the program and they compared favorably 
with known previous results. In the case of pure forced convec- 
tion, laminar flow with a negligible volume heat source, the ve- 
locity profile and the Nusselt number checked exactly with Dr. 
Ojalvo’s results for the same problem with no volume heat sources. 
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The results for turbulent heat transfer in combined forced 
and free convection with volume heat sources showed that an in- 
crease in Prandtl number had a smaller effect than an increase 
in Re* on increasing Nusselt number. The volume heat source 
parameter P had negligible effect on Hu „ Rayleigh number had 
no effect on Nu. The pressure-drop parameter C increased approx- 
imately an order of magnitude as Re* increased an order of mag- 
nitude. Prandtl number Pr, Ra, and P had no effect on C for 
Ra <C 625. When Ra > 625 and P > 1 , increasing Pr decreased C 
and increasing Ra increased C; and for Ra > 625 and P < 1 , in- 
creasing Pr increased C and increasing Ra decreased C. Decreas- 
ing P for Ra > 625 decreased C, Increasing Pr had less effect 
than increasing Re* on decreasing the mixed -mean- to- wall temper- 
ature difference cj> Vl . Increasing Ra had negligible effect on 
• Decreasing P increased cj) fn . 

The velocity in the center of the tube lowered and became 
less positive a 3 Pr, Re*, or volume heat sources were Increased. 
When F < 1 , increasing Ra flattened the velocity profile. 

Increasing Re* or Pr flattened the temperature difference 
profile (|» . Increasing Ra had negligible effect on "(f) . De- 

creasing P increased (f> and when F - 1 , the temperature differ- 
ence profile passed through zero. 
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CHAPTER I 



INTRODUCTION 



The first theoretical investigations into combined forced 



and free laminar flow in a vertical tube were conducted by 
Os tr oumov 10 "' and Hallman^. In addition to predicting velocity 
and temperature profiles. Hallman extended the analysis to in- 
clude calculation of Nusselt numbers and pressure drops. He 
considered cases in which volume heat sources were either pres- 
ent or not present. He also conducted an experimental investi- 
gations ..of the problem for laminar flow. 



The first theoretical investigation into combined forced 
and free turbulent convection in a. vertical tube was conducted 
by Ojalvo and Grosh^. The analysis did not consider cases with 
volume heat sources present. All cases of laminar heat trans- 
fer, pure forced convection or combined forced and free convec- 
tion, checked exactly with known results. Cases of pure forced 
convection turbulent flow were not in agreement with known re- 



sults apparently because too high a value of eddy viscosity was 
used in the buffer zone and in part of the turbulent core. The 



^'Superscript numbers 
the Bibliography. 



•efer to similarly numbered 



references 



in 
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results of combined forced and free turbulent convection were 
presented. 

Sackett-*-3 presented data for combined forced and free tur- 
bulent convection in a vertical tube with a new relationship for 
the eddy diffusivity of momentum. The results for turbulent 
flow cases indicated some improvement but did not compare favor- 
ably with Dr. OjaJ.vo ' s results. The most likely cause was attri- 
buted to the reworked computer program. 

Jackson^ conducted an investigation to improve the relation- 
ship for the eddy diffusivity of momentum. The results of his 
study are used in the present study. 

The present analysis investigates the transfer of uniform 
thermal energy from the wall of a round tube and. uniform ther- 
mal energy from volume heat sources to a fluid flowing vertically 
upwards in the tube. Both forced and free convection exist in 
the steady turbulent flow. Only the case of fully developed flow 
is considered and compared with known results. 

In reworking the original computer program, an error was 
found that would affect Dr. Ojalvo’s results for combined forced 
and free turbulent convection. This error did not occur in 
Sackett’s program, thus the apparent cause for the differences 
in their results. 
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'j?he analysis has practical applications in the fields of 



nuclear reactors and heat exchangers. Its 
the effect of uniform volume heat sources 
and free turbulent convection problems. 



3 solution will show 
on combined forced 
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CHAPTER II 
ANALYSIS 



The problem to be analyzed if? combined forced and free 
turbulent convection with uniform volume beat sources in a 
vertical circular tube whose axis is parallel to the direction 
of the body force. There is to be a net through-flow. 



Assum p tion s 

In addition to the description of the problem given above., 
the following assumptions are made: 

1. Axial symmetry exists for the momentum and heat 
transfer. 

2. All fluid properties are constant , except density 
in the expression for body force. A mean density 

is used for all other density terms. 

3. Viscous dissipation and axial heat conduction are 
negligible compared with the heat conduction in 
the radial direction. 

4. There is uniform wall heat flux. 

5. The velocity and temperature profiles are fully de- 
veloped. There are no radial or angular velocity 
components . 

6. There is single-phase flow. 



"The Nomenclature is given In Appendix A. 
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7. Flow Is steady, turbulent and Incompressible. 

8. The eddy dlffusivity of momentum is given by Jack- 
son's equations. (See Appendix B, ) 



Basic Equa tlona 

The basic equations employed in this analysis are the 
continuity, Navier-Stokes, and energy equations in cylindrical 
coordinates^. An equation of state is also used. 

On the basis of the above assumptions and description of 
the problem, the continuity equation 

v • ( p V ) = 0 
hr ^ 



reduces to 

= o , 

^X y 

the Navier-Stokes (Momentum) equations 

^ fr ~ ® * /• + ~3, y j ? ( v f v ) 

reduces to 




and the energy equation 



= v - Icvt + f 

reduces to 




q' U — -L AV r T !• • n C-- c.;l bt \ 

l r hr l > vl d», * H ]j f : ) 



(i) 



(2) 

(3) 



(4) 
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The equation of state is developed from an expansion of ^ in 
a Taylor series in t about the reference temperature t w , and 

is : 8 

? = ew[i-p(t-tw)] • (5) 

D evelopment of Equati on s 

Equation (l) indicates that 

u = u (r), ( 6 ) 

and equations (3) indicate that 

P = P (x) . ( 7 ) 

If equations (5), ( 6 ), and (7) are utilised, equation (2) can 
be written as 




The boundary condition of uniform wall heat flux plus 
the assumptions of constant specific heat, uniform volume heat 
sources, and a fully developed temperature profile require that 



it _ 



R 



(a constant). 



( 9 ) 



Equation (9) is developed in Appendix C. A new variable 0 
is introduced and defined as 



O = 0(r) = t(x,r) - 1 (X 





(10) 
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In terms of equation (10) , equations (4) and (8) become. 



respectively. 



f^ ufl “ f= T 



( 11 ) 



and 



-L £_( r Vu fiaial du ( Y \ , p A £ i)(r) - dp(x.)j... { o 1 
r A. [ 1 r Q j r [r '1 \wV cj r 'J' fl ~ d ' x i \w q 



r <H L 5c ?c 



' w ' & 



c/x 



c/C 



(12) 



Although the density ^ in equation (12) is actually a 
function of x, by considering the solutions of equations (11) 
and (12) to obtain u and Q at a fixed value of x, we can con- 
sider constant®. Our assumption of fully developed flow 
and heat transfer is consistent with these results. Each side 
of equation (12) may be set equal to some constant since vie 
have separated the variables^. Thus, 



and 



d? , o A - 
dx Vw % 



_ 3 1 u mju C (a constant ) 

D? 3c 



(13) 



f-i( r [i + ¥ld?) + f w? i 0 = - a y^ t£L • w 

The pressure-drop parameter C in these equations takes on the 
value of unity for the special case of pure forced-convection 
laminar low, as described by Hallman3 . 
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In order to solve the set of problems described, equations 



(ll) and (l4) are nondimensiona 
dimensionless terms: 

)|e 2 r/ D 

(j) s i$k e/i" d ? 

U — f,-i 

F = f m u m Cy> 1 /^"' 

Equation (19) is discussed in A 
Thus, equation (ll) becomes 



..ized by using the following 



radius 


(15) 


t empe ra ture di f f er enc e 


( 16 ) 


velocity 


(17) 


Rayleigh number 


( 18 ) 


volume heat source 


(19) 


pendix F. 





_L A 

y\ ^ 




- in ] .di \ 

^ m. d 7 ] / 



and equation (l4) becomes 



- f FU ~ f 



( 20 ) 



r 1- (y\ 


1 f 


6m 


41/ 


\ _ Ra <t> 


- sc : . 


I) d7|l 1 




7 / _ 


c/ 


J ” fp 


In accordance 


with as 


sumption (8), 


we will use 


_isd “ g" 




6 


05 


_ r / 




“ 




Tr 2 


| 


a/ -- e x y> lot is! 2 / 


?~p~) 



( 21 ) 



( 22 ) 



as given by L-ykoudJ s 
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Equation (20) may thus be written as 



_L A 

3^ 




7/ 



d 4> \ - 

00 y 



F U - l , L 



( 23 ) 



n 



Bo und s ry Condlt Ion s 

The following boundary conditions are used: 



uo) - o m 

(J> 0) — 0 (25) 

dJlM = u'( 0 ) - 0 (26) 

d >1 

cL^.Xo) ^ (j) / ( o) “ 0 (27) 

di) 



These boundary conditions come from the physical problem. 
Equations (24) and (25) state that the fluid velocity and tem- 
perature at the v/all of the tube (Tj z l) are equal, respec- 
tively, to the wall velocity and wall temperature. Equations 
(26) and (27) come from the fact that heat and momentum are 
not transferred across the center line (yj = Q) of the tube, 
due to axial symmetry, resulting in a zero slope for the temper- 
ature and velocity profiles at the center line7. 

Dis cussion o f Equa tions 

Equation (21), the momentum equation, and equation (23), 
the energy equation, are to be solved by a 1620 IBM digital 
computer for U and ([) as functions of 7| • All boundary con- 
ditions are to be satisfied. The value of the pressure-drop 
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parameter C may also be obtained if we use the following form 



of the continuity equation: 



Urn ~ d U 7) dTj 



(28) 



Three dimensionless parameters in our equations are the 
Prandtl number Pr, the Rayleigh number Ra , and the volume heat 
source parameter F. The Prandtl number determines the rela- 
tionship between the velocity and temperature distributions. 
For the special case of Pr ~ I, both profiles are the same. 
Rayleigh number Ra - Gr Pr where the Grashof number Gr repre- 
sents a ratio of buoyant forces to viscous forces 12 . Thus, 
free-convection predominates for large Rayleigh numbers and 
for the special case of Ra - 0, pure forced-convection is rep- 
resented. The significance of variations in the magnitude of 
the volume heat source parameter F is discussed in Appendix F. 

Empirical equations for are given in Appendix B. 



They are : ^ 






c /?eV5' oo f or q.< -L , ( 3 i) 

D ^ 

where Ke z 0 U . , friction Reynolds number 
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Thus, a fourth parameter, the friction Reynolds number Re*, 
is introduced. In this analysis, when Re* - 0, the ratio 
is set equal to zero, thus assisting in checking the validity 
of the program by comparison with laminar flow problem results. 

Extensio n of Res u lts 

The solution of equations (21), (23), and (28) will yield 
U and ([) as functions of 7| for a given set of values of the 
parameters: Ra, FT, Re*, and F. The pressure-drop parameter 

C will also be obtained. 

These results may be extended by calculating the following 
useful quantities: the dimensionless mixed-mean-to-wall tem- 
perature difference ; the Nusselt number Nu; and the Rey- 

nolds number Re. 

Hallman3 gives equations for the first two of these quan- 
tities. They are: 

K - ^ (34) 

and 

A/m E -£■ (l -F) ; (35) 

see Appendix D. 

An expression for the Reynolds number is derived In 



1.4 



Appendix E. 


The result is 

+ / p r \2 


Re = 


- f K e ) 

SC+RaQm/^F ' 


where the -{- 


sign is for upward floitf and the - sign is for 


net downward 


flow. 
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CHAPTER III 



METHOD OF SOLUTION 



The method of solution is based on the use of a digital 
computer. This method will decrease considerably the amount 
of hand calculations required in the solution of our problems. 
The possibility of errors are thus reduced. This method of 
solution will also allow us to investigate many different 
problems and compare some of these results with known analy- 
tical and experimental data. 

The independent parameters used are Ra, Re*, Pu, and F. 

V 

The Rayleigh number is chosen to measure the extent of the 
free-convective effect on the flow. The friction Reynolds 
number determines the value of -~p- as a function of Ij . 

The Prandtl number determines the value of <T~ . The volume 
heat source parameter is chosen to describe the thermal ener- 
gy convected downstream to the heat generated in the fluid. 
With these input quantities fixed, the method of solution pro- 
ceeds as follows: 



is assumed as an initial guess for the velocity profile. This 



1 . 




( 37 ) 



equation satisfies the continuity equation 
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and the boundary conditions U(l) ~ 0 S and U ! (Q) ~ 0, It is 
the limiting case of the parabolic profile in pure forced- 
convection laminar flow. 

2. Equation (23) is Integrated with the aid of boundary 



condition (ff (Q)~ 0 to obtain 

' ■'* _ 



r = pf/q U:'n<b \ _ <L 2 l 

^ ">1 fl -!- <r-Pr taA fl +<T? 



(38) 



^ [| -h <r?r [| + <T p r laj 

3. Equation (38) is integrated using boundary condition 

<j> a) = 0 to obtain 

as a function of 3^ . 

b . Equation (21) is Integrated using boundary condition 



(j) = p 



(39) 



U* (0) z 0 to obtain 

l^D.-lol$2\.A2L 



Ui 



2 



v-j:uC_ 



[h- 






[l f %r] 



(‘to) 



The pressure-drop parameter C must have an output numer- 
ical value for the computer. For the limiting case of pure 
forced-convection laminar flow C - 1. 

5. Equation (40) Is integrated subject to boundary condi- 
tion U(l) ~ 0 to obtain 

V, --pv'6^ {41) 

as a function of "y^ . 

6. The continuity equation (28) is checked to see if 
equation (4l) satisfies it. If U2 does not satisfy this couat: 



ion. 
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the value of C is changed until it is satisfied. 

7. The function Ug is checked against Uj by comparing 
U' 2 (l) with U'j^l) to see if they are within 0. 1$ of each other. 

8. Equations 



are used to calculate (j) M , Nu and Reg respectively. If the 
check in Step 7 is satisfactory. The results are then printed 
out as given in Step 10 below. 

9. If the check in Step 7 is not satisfactory, U2 is used 
as an initial assumption, replacing Ui in going through the 
procedure again, starting at Step 2. 

10. Print out 



2 1 

A >u- * 




) 



and Re 



Title 



Pr obi era Numb e r 



F 



Ra 



Re 



G 



d 




1 



V 

0 



V 

0 



0 



c 




Re 
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VI 2 and are both printed out to ensure that the criterion 

for checking them in Step 7 is valid. 

These steps may be shown in a simplified flow diagram for 
the calculation. Figure 1. 
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CHAPTER IV 



COMPARISON OP RESULTS 



Table 1 lists the following eal ciilated results: C, <$) m . 

Re, and Nu for all problems solved by the 1620 digital computer. 



TABLE 3. 









RESULTS FOR 


, ALL COMPLETED PROBLEMS 








Inpu t Pa r ame t er s 


Calculated 


Results 


Problem 














No. 


Ra 


5j{ 

vD 

Cm 


Pr 


F 


6 


Re 


Nu 








A. Laminar 


flow, pure forced convection 


403B 


0 


0 


--- 


103 


1.0 - 915.85 




4.36 








B. Laminar 


flow, combined fore 


ed and free 


402B 


81 


0 


— 


10 


2.4 - 9.62 


-- 


3.74 








C. Turbulent flow, pure forced 


convection 


379 


0 


103 


1 


.5 


7.3 .036 


17000 


55.4 


380 


0 


10 3 


1 


1.0 


7-3 -.0034 


17000 


0.0 


383 


0 


103 


1 


10.0 


7.3 -.71 


17000 


50.3 


385 


0 


103 


10 


.5 


7.3 .01 


17000 


196.0 


386 


0 


103 


10 


1.0 


7.3 -.00053 


17000 


0.0 


389 


0 


103 


10 


10.0 


7.3 -.194 


17000 


185.0 


391 


0 


10 3 


100 


.5 


7.3 .002 


17000 


983.0 


392 


0 


10 3 


100 


1.0 


7.3 -.000064 


17000 


0.0 


395 


0 


103 


100 


10.0 


7.3 -.038 


17000 


952.0 


397 


0 


10^ 


1 


.5 


53.9 .005 


232000 


400.0 


398 


0 


10 £ ' 


1 


1.0 


53.9 -.00017 


232000 


0.0 


400B 


0 


300 


1 


10.0 


2.8 - 2.07 


4030 


17.4 


4oi 


0 


10^ 


1 


10.0 


53.9 -.093 


232000 


385.0 


403 


0 


10^ 


10 


.5 


53.9 .0011 


232000 


1750.0 


4 04 


0 


10? 


10 


1.0 


53.9 -.000019 


232000 


0.0 


407 


0 


10^ 


10 


10.0 


53.9 -.02 


232000 


1720.0 


409 


0 


10 4 


100 


.5 


53.9 .00025 


232000 


7900.0 
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Table 1 (Continued) 



Problem 


Input 


Parameters 




Cal cu.lated 


Results 


No. 


Ra Re* 


Pr 


p 


c 




Re 


Nu 




C. 


Turbulent flow,, pure forced convection (Continued) 


4l0 


0 


10 


100 


1.0 


53.9 


-.000002 


232000 


0.0 


413 


0 


10 il 


100 


10.0 


53 . 9 


-.0046 


23 2000 


7850 . 0 


415 


0 


10 5 


1 


.5 


430.0 


.00063 


2900000 


3150.0 


4l6 


0 


io5 


1 


1.0 


430.0 


-.000021 


2900000 


0.0 


419 


0 


105 


1 


10.0 


430.0 


-.012 


2900000 


3040.0 


421 


0 


lo5 


10 


.5 


430.0 


.00017 


2900000 


11700.0 


422 


0 


io5 


10 


1.0 


430.0 


-.0000034 


2900000 


0.0 


425 


0 


105 


10 


10.0 


430.0 


-.0031 


290^000 


11500.0 


427 


0 


io5 


100 


.5 


430.0 


.0001 


2900000 


18600.0 


428 


0 


10-5 


ICO 


1.0 


430.0 


.-.0000014 


2900000 


0.0 


429B 


0 


10 4 1 103 

I). Turbulent flow. 


53.9 -'10.34 
, combined forced 


232000 
and free 


386.0 


55 


81 


10 3 


1 


.1 


5.8 


.067 


1.6600 


53.1 


56 


81 


10 3 


1 


.5 


7.1 


.036 


17000 


55.4 


57 


81 


10 3 


1 


1.0 


7.3 


“.0036 


17000 


0.0 


60 


81 


10 3 


1 


10.0 


7.5 


-.718 


17000 


50.0 


61 


81 


10 3 

10 3 


10 


.1 


6.9 


.018 


16900 


191.0 


62 


81 


10 


.5 


7.3 


.01 


17000 


196.0 


63 


81 


10 3 


10 


1.0 


7.3 


-.00054 


1.7000 


0.0 


66 


81 


10 3 


10 


10.0 


7.4 


-.194 


17000 


185.0 


6 7 


81 


109 

10 3 


100 


.1 


7.2 


.0037 


17000 


970.0 


68 


81 


100 


.5 


7.3 


.002 


I.70OO 


983.0 


69 


81 


10 3 

10? 


100 


1.0 


7.3 


-.000064 


17000 


0.0 


72 


81 


100 


10.0 


7.3 


-.037 


17000 


952.0 


73 


81 


10 


1 


.1 


53.7 


. 0091 


232000 


394.0 


74 


81 


10 


1 


.5 


53.9 


.005 


232000 


400.0 


75 


81 


10 4 

10 ll 


1 


1.0 


53.9 


-.00017 


232000 


0.0 


78 


81 


1 


10.0 


53.9 


-.093 


232000 


385.0 


79 


81 


10 

10? 


10 


.1 


53.8 


.002 


232000 


1740.0 


80 


81 


10 


.5 


53.9 


.001 


232000 


1750.0 


81 


81 


10, 

10? 

1°; 

10, 


10 


1.0 


53.9 


-.00002 


232000 


0.0 


84 


81 


10 


10.0 


53.9 


-.02.1 


232000 


1720.0 


85 


81 


100 


.1 


53.9 


.00045 


232000 


7890 . 0 


86 


81 


1.00 


.5 


53.9 


.00025 


232000 


7900.0 


87 


81 


10/ 

io ; > 


100 


1.0 


53.9 


- .000002 


232000 


0.0 


90 


81 


100 


10.0 


53.9 


-.0045 


232000 


7850.0 
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Tabic 1 (Continued) 

Input Parameters 

Problem 



Calculated Results 



No. Ra 


Re* 


Pr 


p 


c 


Qrn 


Re 


Nu 




D. 


Turbulent flow, combined 


forced and 


free (Continued) 


91 


81 


10? 


1 


.1 


427.0 


.0011 


2920000 


3120.0 


92 


81 


10 = 


1 


.5 


430.0 


.00063 


2900000 


3150.0 


93 


81 


10 5 


1 


1.0 


430.0 


. 000021 


2900000 


0.0 


96 


81 


10 5 

r~ 


1 


10.0 


430.0 


-.012 


2900000 


3040.0 


97 


81 


10 ? > 


10 


.1 


427.0 


. 0003 


2920000 


11700.0 


98 


81 


10 J 


10 


.5 


430.0 


.00017 


2900000 


11700.0 


99 


81 


3.0 5 


10 


1.0 


430.0 


-.0000034 


2900000 


0.0 


102 . 


81 


105 


10 


10.0 


430.0 


-.0031 


2000000 


11500.0 


103 


81 


106 


100 


.1 


427.0 


.00019 


2920000 


18600.0 


104 


81 


lo5 


100 


.5 


430.0 


.0001 


2900000 


18600.0 


105 


81 


io5 


100 


1.0 


430.0 


-.0000014 


2900000 


0.0 


108 


81 


lo5 


100 


10.0 


430.0 


-.0019 


2900000 


18400.0 


164 


625 


iq3 


1 


.5 


6.0 


.036 


16700 


55-4 


165 


625 


10 3 


1 


1.0 


7.4 


-.0035 


17000 


0.0 


169 


625 


10 3 


10 


.1 


4.0 


.019 


16200 


191.0 


170 


625 


10 3 


10 


.5 


7.0 


.010 


17300 


196.0 


171 


625 


10 3 


10 


1.0 


7.3 


-.00056 


17000 


0.0 


174 


625 


10 3 


10 


10.0 


7.6 


-.194 


32000 


185.0 


175 


625 


10^ 


100 


.1 


6.7 


.0037 


16900 


970.0 


176 


625 


10-^ 


100 


.5 


7.2 


.002 


17100 


983.0 


177 


625 


10 3 


100 


1.0 


7.3 


-.000064 


17000 


0.0 


180 


625 


10 3 


100 


10.0 


7.4 


-.038 


18800 


952.0 


181 


625 


10 


1 


.1 


52.2 


.0091 


231000 


394.0 


182 


625 


10 


1 


.5 


53.7 


.005 


232000 


4oo.o 


183 


625 


10 


1 


1.0 


53.9 


- . 00017 


232000 


0.0 


187 


625 


10 l 


10 


.1 


53.5 


.002 


232000 


1740.0 


188 


625 


10 4 


10 


.5 


53.8 


.0011 


232000 


1750.0 


189 


625 


10 4 


10 


1.0 


53.9 


-.000019 


232000 


0.0 


192 


625 


10 4 


10 


10.0 


53.9 


— .021 


233000 


1720.0 


193 


625 


10 4 


100 


.1 


53.8 


.00046 


232000 


7890.0 


194 


625 


10 4 


100 


.5 


53.9 


.00025 


232000 


7910.0 


195 


625 


10 4 


100 


1.0 


53.9 


-.000002 


232000 


0.0 


198 


625 


10 4 


100 


10.0 


53.9 


-.0046 


232000 


7850 . 0 


199 


625 


10 5 


1 


.1 


427.0 


.0011 


2920000 


3120.0 


200 


625 


io5 


1 


.5 


430.0 


.00063 


2900000 


3150.0 


201 


625 


105 


1 


1.0 


430.0 


-.000021 


2900000 


0.0 


205 


625 


10 5 


10 


.1 


427.0 


.0003 


2920000 


11700.0 
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Table 1 (Continued) 

Input Parameters Calculated Results 

Problem 



No. 


Ra 


Re* 


Pr 


p 


c 


YW 


Re 


Nu 




D. Turbulent flow., combined 


forced and 


free (Continued) 


206 


625 


10* 


10 


.5 


430.0 


.00017 


2910000 


11700.0 


207 


625 


10 5 


10 


1.0 


430.0 


-.0000034 


2910000 


0.0 


210 


625 


io5 


10 


10.0 


430.0 


-.0031 


29-10000 


11500.0 


211 


625 


10 5 


100 


.1 


427 .0 


.00019 


2920000 


18600.0 


212 


625 


10 5 


100 


.5 


430 . 0 


.0001 


2910000 


18600.0 


213 


625 


io5 


100 


1.0 


430.0 


-.0000014 


2910000 


0.0 


216 


625 


io 5 


100 


10.0 


430.0 


-.0019 


2910000 


18300.0 


327 


10J 


103 


1 


1.0 


8.3 


-.0035 


3 7400 


0.0 


330 


10 i 


103 


1 


10.0 


27.6 


-.73 


27600 


48.7 


332 


10 4 


103 


10 


.5 


1.6 


.01 


15700 


196.0 


333 


10 1 


IO 15 


10 


1.0 


7.5 


-.00054 


17100 


0.0 


336 


10 


10 3 


10 


10.0 


12.7 


-.194 


18700 


I85.O 


338 


10 2 * 


io3 


100 


.5 


6.2 


.002 


167OO 


984.0 


339 


i°2 


Io3 


100 


1.0 


7.4 


-.000072 


17000 


0.0 


342 


10 £ 


IO 3 


100 


10.0 


8.4 


-.038 


17300 


948.0 


343 


10 £ 


io4 


•1 


.1 


27.4 


.0091 


223000 


393.0 


344- 


10? 


10 t 


1 


.5 


51.1 


.005 


230000 


400.0 


345 


10 a 


10 4 


1 


1.0 


54.0 


-.0002 


231000 


0.0 


348 


io J 


10 a 


1 


10.0 


56.7 


-.094 


232000 


383.0 


349 


io , 


10* 


10 


.1 


47.8 


.002 


230000 


1740.0 


350 


10 4 


107 

jl 


10 


.5 


53.2 


.0011 


232000 


1750.0 


351 




10* 


10 


1.0 


53.9 


-.000019 


232000 


0.0 


354 


10 h 


10 4 


10 


10.0 


54.5 


-.021 


232000 


1720.0 


355 


10 a 


10* 


100 


.1 


52.5 


.00045 


231000 


789O.O 


356 


10 ? 


10 


100 


.5 


53.7 


.00025 


232000 


7910.0 


357 


10 1 


1 °5 


100 


1.0 


53.9 


-.000002 


232000 


0.0 


360 


10* 


10 4 


100 


10.5 


54.0 


-.0046 


232000 


7850.0 


36l 


10 4 

I . 


10 5 


1 


.1 


424.0 


.0011 


2920000 


3120.0 


362 


10 4 


105 


1 


.5 


430.0 


.U0063 


2900000 


3150.0 


363 


10 J 


105 


1 


1.0 


430 . 0 


-.00006 


2900000 


0.0 


3 66 


lof 


105 


1 


10.0 


430.0 


-.012 


2910000 


3050 . 0 


367 


10! 


io5 


10 


.1 


426.0 


.00003 


2920000 


11700.0 


368 


10, 


lo5 


10 


r— 

0 


430 . 0 


.00017 


2910000 


11700.0 


369 


1°, 


io5 


10 


1.0 


430.0 


-.0000034 


2910000 


0.0 


372 


1°,; 


105 


10 


10.0 


430 . 0 


-.0031 


2910000 


11500.0 


373 


10 


IO 5 

r— 


100 


.1 


427.0 


.00019 


2920000 


18600.0 


374 


i° 


10 


100 


.5 


430.0 


.0001 


2910000 


18600.0 


375 


i° 


10 n 


100 


1.0 


430.0 


-.0000014 


2910000 


0.0 


378 


10 4 


IO 5 


100 


10.0 


430.0 


-.0019 


2910000 


18400.0 



Part A consists of the problem in laminar flow* pure forced 
convection with negligible volume heat sources. Part B con- 
sists of a problem in laminar flow, combined forced and free 
convection with volume heat sources. Part C consists of 28 
problems in turbulent flow, pure forced convection with volume 
heat sources. Fart 1) consists of 100 problems in turbulent 
flow , combined forced and free convection with volume heat 
sources . 

The IBM solutions to all problems are an attachment to 
the present analysis. Copies of the Fortran program, the data 
cards, and the solutions to all problems on IBM cards are in 
the files of Dr. Ojalvo, 

Discus si o n 

The velocity profile and the Nusselt number of the laminar 
flow problem in part A, Table 1, checked exactly with the re- 
sults of Ojalvo® and Hallman®. A negligible volume heat source, 
F=1000, was used in this problem to check the validity of the 
program. The laminar flow problem, combined forced and free 
convection with volume heat sources of part B, Table 1, com- 
pared favorably with Hallman. It also compared favorably with 
Ojalvo even though no volume heat sources were present in his 
analysis . 
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Problem 397 was chosen as typical of the problems in part 
C. This problem was compared with test problem t'29B and the 
universal velocity and temperature profiles of Eckert and 

p 

Drake . 



The values of U and (j) m in problem 429B with negligible 
volume heat sources, FnlOOO, were checked by using the defini- 



tions' 



U 



y 



yhwTpyX/ , 



and 



(At) 



+ 



t - t 



w 



(42) 



(43) 



(44) 



with the definitions of the various dimensionless quantities 



to derive the comparison equations 

U + = -f g, (17) 

K e 



(45) 



and 



y + - .!. Rz ( i -~n) , 

(zt + )~ 1 jPrlsL ((f)) 

” ‘I (l-F) W 



(46) 

(47) 



Figure 2 shows the comparison of u ^ and - (ot" ! ) for pro:- 
lem 429B with the universal velocity and temperature profiles. 
The values of velocity coincide almost exactly up to y - *' =3 00, 
then they start to diverge slowly below the universal profile. 
The values of temperature fall slightly above the universal 
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temperature profile. The •values ol 
with those of problem 4-29B except 



problem 397 concide exactly 
the temperature starts to 



diverge above the values 



of problem 429B at y r ^300. 



Volume 



heat sources had negligible effect on the velocity and temper- 
ature profiles. 



The differences in the 



problem 3 of Dr. 



Oji 



, 8 . 

live are 



values of these two problems and 
attributed to the improved expres- 



sion for the relative viscosity. 



A decrease in F from 1000 to 



that is an increase in 



the 



volume heat sources ^ had very 



little effect on the profiles 



All of the problems in part C with a Prandtl number of 1 



or 10 and with a volume heat source parameter F-10 are plotted 
in Figure 3, which shows the variation of Nusselt number with 
Reynolds number. The Dittus-BocO tea* equation? for pure forced- 
convection turbulent heat transfer is also plotted. In both 
Prtmdtl number comparisons the results of the present analysis 



compare almost exactly with the c 
equation. Volume heat sources ha 
variation of Nusselt number with 
Figures 4 through 10 shun r-e 
bulent flqvg combined fore;-.;'' and 
number is plotted against Ray lei 



urves of the Dittus-Boelter 
d nog. 1 igible effect on the 
Reyno 1 d s numb e v , 
cults for the problems in tun- 
free convection. The Nusselt 
h number in Figure 4. The 



n f y, 



Nun soil; number is increased less than an order of magnitude as 
Re* is increased an order of magnitude. Increasing the Prandtl 
number has less effect on increasing the Nussult number them 
increasing Re*. Increasing the Rayleigh number had no effect 
on Nusselt number. Decreasing F had a very little increasing 
effect on Nu. 

Several test problems were run to determine the limitations 
of the computer for Re*=10^. No solutions were obtained on the 
IBM 1620 computer for Re* =10 3., F=.l or less., and Ra-~625 or more. 
Only by decreasing Ra or increasing either Re* or P 'would the 
computer give a solution. A large Ra in combination with a 
small F (a large volume heat source) would physically appear 
to be opposing each other. Convective forces predominate in 
a large Ra., whereas with a small F, considerably more heat is 
being generated than is being converted downstream. 

Figure 5 shows the pressure-drop parameter C plotted 
against Rayleigh number. C increased approximately an order of 
magnitude as Re* increased an order of magnitude. For Rayleigh 
number less than 625', Ra* F, and Pr had no effect on C. For 
values of Ra over 625 and F'tl, increasing Pr decreases C and 
increasing Ra increases C. If F<1, increasing Pr increases C 
and increasing Ra decreases C. These curves show that C changes 
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slope for a value of P approximately equal to 1. For a given 
Re* and Pr the point at which C changes slope is a function of 
Ra. For small values of P and Re*, and large values of Ra, C 
would appear to go negative verifying Hallman's results for 
laminar f lovP „ As Re* is increased the pumping action of large 
heat sources is reduced considerably. 



Figures 6, 7, and 8 show the variation of mixed-mean- to- 
wall temperature difference with Rayleigh number. Increasing 
Re* an order of magnitude approximately decreases (£) an order 
of magnitude. Increasing Ra has no effect on (j) . Increasing 
Pr has less effect on decreasing fj) m than increasing Re*. 
Finally increasing F has the most effect on decreasing <j) . 

As shown in Figure 8, (j) passes through sere for F approxi- 
mately equal to 1 . If (j) m is plotted against l/F, as is shewn 
in Figure Sh it is seen that a variation in Ra, Pr, or Re* has 
little effect on this crossover point. Hallman's curve for 
laminar flow, Ra=0 is plotted for comparison. 

When (j) ]n -0, Nu=oo, according to equation D~3. Figure 10 
shows Nusselt number plotted against l/p for Ra-0, Re*=lo3, and 
Pr of 1 and 10. An infinite discontinuity occurs whenever 
(j) m =0. From equation D-3 it is seen that Nu=0 for F=l. Hall- 
man's plot for laminar flow, Ra=0 is shown for comparison. 
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Representative velocity and temperature profiles are shown 

' V 

in Figures 11 through 21. 

The effect of varying Re* is shown in Figures 11 through 
1 6 for Pr~-10, Ra=625> and for various values of F. As Re* 

O [{ 

increases from 10-’ to 10 , the velocity in the center of the 
tube becomes less positive and the velocity near the wall be- 
comes more positive. As Re* is increased, the temperature 
profiles flatten out and the slope becomes higher at the wall. 
Decreasing F, as shown in Figure 21, flattens the velocity pro- 
file, whereas the temperature profile decreased considerably 
and passes through zero at F=l. The effect of increasing F on 
the velocity and temperature profiles was not as is predicted 
in appendix F. The action is due to the increased values of 
Re* investigated in the turbulent flow problems. 

Figures 17 and 18 show the effect of various Prandtl nurn- 

o h 

bers for Re^rlO-*, RazlO , and F=10. As Prandtl number increases, 
the velocity and temperature profiles decrease more in the cen- 
ter of the tube than near the wall, thus becoming more flat. 

The effect of a variation of Ra for Re*=10^, Pr=lO, and 
P’=,5 is shown in Figures 19 and 20. As Ra is increased from 0 
to 10 4 there is very little decrease in the velocity and tem- 
perature profiles at the center of the tube. 
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CHAPTER V 



SUMMARY AMD CONCLUSIONS 

The effect of the parameters Re*, Pr, Ra, and F on Nu, C, 
$ nw U, and (?) profiles is summarized in Table 2. 



TABLE 2 

SUMMARY OP EFFECTS IN COMBINED 
FORCED AND FREE CONVECTION 

Change on Nu C ([Vn u $ 

for , Ra<625 Ra>625 Velocity Temperature 

profile difference 
profile 

F> 1 

increasing 



Re* 


increase 


increase 


increase 


decrease 


flattens 


flattens 














curve 


curve and 
makes (j) 
approach 
















zero 


Pr 


increase 


no 


effect 


decrease 


decrease 


lowers 


same as 














velocity 
in cen- 
ter of 
tube , mak- 
ing it 
less pos- 
itive 


Re* 


Ra 


no effect 


no 


effect 


Increase 


slight 


increase 


decreases 












decrease 


in center 


in center 












at high 
Ra 


of tube 


of tube 
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Table 2 (Continued) 



Change on Nu 

for 

F=1 

increasing 


C 

Ra<625 


Ra>625 




U 

Profile 


4> 

Profile 


Re* 


0 


increase 


increase 


decrease 


flattens 

curve 


negligible 

effect 


Pr 


0 


no effect 


decrease 


decrease 


no effect 


negligible 

effect 


Ra 


0 


no effect 


no effect 


no effect 


no effect 


no effect 



F< 1 

increasing 



Re* increase increase increase decrease flattens same as 

curve F > 1 



Pr increase no effect increase decrease flattens same as 

curve F > 1 - 



Ra no effect no effect 



decrease no effect flattens 

curve 

slightly 



negligible 

effect 



decreasing 

F slight no effect decrease increase flattens increase, 
increase curve passing 

through 
zero, F=1 



Conclusions 



1. All test problems checked sufficiently close to previous 
known results to insure the validity of the program. 

2. The new expression used for the relative viscosity enabled 
the universal velocity and temperature profiles to be approached 
almost exactly in pure forced-convection turbulent heat transfer 

y problems with negligible volume heat sources. 

3. Volume heat sources had negligible effect on the velocity 
profile, but had considerable effect on the temperature profile 
as would be expected from the development of the equations. 

4. The program was unable to solve problems involving high 
Rayleigh numbers and volume heat sources for low values of Re*. 

5. Increasing the various input parameters tend to flatten 
the velocity and temperature profiles as would be expected. 

6. The effect of volume heat sources was reduced considerably 
by the influence of turbulent flow. 
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Figure 1 . Simplified Flow Diagram of the Calculations 
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Nu - 0.023 (Re ) (Pr) 

Results of Present analysis 
for P 10 



4 
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NUSSELT NUMBER, Nu 
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Figure 4. Variation of Nusselt Number 
with Rayleigh Number 
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PRESSURE DROP PARAMETE 
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Velocity Profiles, F=l, 
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Figure 13. 
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Figure 14. Temperature Difference Profile,, 
Frl j, PrrlOj Rar625 
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Figure 15- Velocity Profiles, F - 10, 
Pr - 10, Ra - 625 
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Figure 16. Temperature Difference Frofile 
*‘=10, Pr=10, Ra=625 
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Figure 17. Velocity Profiles., FrlO, 
Re--] o3 j Pie -.10 4 
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Figure 18. Temperature Difference Profiles, 
F-10, Re*=10 3 , Ra=10 2 ' 
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Figure 19. Velocity Profiles, F=5, Pr=10, 
Pie -r 10 3 
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Figure 20. Temperature Difference Profile, 
F=.5, Pr=10, Re*=10 3 
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APKiHPXa a 



A axial teaipar&turo gradient in fluid f y~ , ry'-ft 
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So 

h 



t 

.> 



C pressure drop parameter* 

d imens i o: il e ,v a '"*■* ' v/ v 



3 c M m ) 

*h specific heat of fluid at constant pressure , BJLJL^ 

^ ' TCPl- 



tube inside diameter 5 ft 

acceleration due to gravity., ft/ sec 2 

dimensional constant in 'Newton's Law, 

4.1? x 10 8 ft lb/ib f hr 2 

convection heat braor f ? r coefficient , 

<J"w/ -0» , E?..,'hr f*; 'V 



k thermal conductivity of fluid, Bu>J/hr ft °F 
p static fluid pressure, lbr/fb- (absolute) 

<j w M heat transfer rate pec sat.? t area at the Fall . 

uiy/hr ft ? 

q nl volttwe hen.t scuz oee., Bve/hr ft 8 
r radial coordinate; I’/’l - r.. ft 

static fluid t 



; -r 



At 


dimer s Lor ’ 1 s-r < 






S' 


tirm, Vr 



1 • i r f 



\ ^ 






fluid velocity parallel to tube axis at radius r, ft/hr 
dimensionless velocity, t Uy^ 
dimens i onl e » s vel oc i t y , u/tt a 



L W 






friction velocity, 
specific volume, ft3/Ib ia 
velocity vector, ft/ hr 

distance measured, along the axis of the tube upwards, ft 
radial coordinate measured from the wall, 1/2 - r., ft 
dimensionless distance from wall, y ^ f - w ' 



t he rma 1 diffus i v ,1 1 y , 



ik = , f**Ar 



t 



I 



eddy diffus ivity of heat transfer, ft 2 /hr 
eddy diffuslvity of momentum transfer, ft^/hr 
dimensionless radius, 2r/D 
radial temperature difference, t - t w , °F 



nr 



d I mens i o nl e a s 



dynamic v5 


SCO 
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of 


,fl 


v i.e. f 1 1 ’j a j 


kinernt J.c 


via 


cos 1 


ty 


of 


fluid, ft 
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o r r 


Ini 


d. 


lb., / ft 3 


ratio of e 
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air 


fun 


ivi 


. » C u 
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wi 


fluid shea 
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tree- 




ln f 


/ p 

/ 1 1< 



dimoTic Lonlfss t,? mweratc re d t f.f erenee, ~ 

•• ' e.n l <ln<7>f)D i 

angular coordinate of fh° cylindrical coordinate system, 
radian': 



SUBSCRIPTS 



H 


heat 


m 


mean 


H 


momenta 


r 


baaed or 


*/ 


radial p 


X 


baaed on 




based on 


t 


differ, u-: 




variable 


! t 


per tun it 


r t f 


per unit 


Gr 


Gras ho f ; 


Nu 


bus celt ! 


Pr 


Prandfcl i 


Ra 


Rayleigh 


Re 


Por y** 1 ^ ; r 

£ t jV i * S -l. .\a 


P L0 W 


P jr ±C o J. 0 1 i 



SUPERSCRIPTS 



NUMBERS 



fe^A^ASjLfil 

" v i . $ 



1 6 y j ~ 



C 7V-/ i\ 



7~VR. 



s r‘.ruoor« 



ji n p^-JD iL 
r r / 6^ /f 

f ^ ra 



Gy Pr 



y 



« 



Ur, 



* i-^V 



58 



APPENDIX B 



THE EDDY DIFFUSIVITY OF MOMENTUM 



An analytical study was made by Jackson^ of the eddy 
diffusivity of momentum in a vertical circular tube. The 
ratio of eddy diffusivity to kinematic viscosity, -~jr— , 

called relative viscosity appeared in a reduced form of the 
Navier-Stokes equation. An empirical fit of experimental data 



was employed to determine the dimensionless velocity gradient 
appearing in the expression for .1,11 

7- 

The result of his study used in the present analysis is 



the following set of expressions for variation of relative 



viscosity with dimensionless radius 



6m - „ VL , 

V - I - .GDC R e * ( I ?//,, -, 0 IS JvJK.iij " I 

in - 9 Re' s'' SO 0 _ I 
V~ 



for 

for 






■f cr 0 £ < _L 

» JO 



The re ;ults compare favorably with the experimental and 
analytical results of others and appears to be an Improvement 
over the Apr expression used in recent studies. These re- 
sults appear as equations (29) j (30), arid (31) in the present 
analysis . 
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APPENDIX 0 



THE LONGITUDINAL TEMPERATDRE GRADIENT 



The heat balance for the steady-state condition in our 
problem states that the rate at which heat is added to the 
fluid is equal to the rate at which heat is eonveeted down- 



stream in the fluid. In symbolic form (see Figure C-l), 



q 1,1 ntf d x _ q " ffD Jx 
y <#• lw 




(c-i) 



where q w " is taken positive if heat is being removed, 
equation (C~l) for 'bty^x gives 



k - r'Phi - %w n 

2 Cy> fo-, X''*ur<)r 

The mean velocity is defined as follows: 

/■Vt. , I 

- i uztrr cl r _ S J urdr 

M " 1 = “ 7VP7T “ 



Solving 



(C-2) 



(C-3) 



from which 




}Ln-> [)} 
8 



(C-4) 



Equation (C-2) can therefore be written as 

K - i-- o Z t j yj. . (c-5) 

Since all quantities on the right-hand side of equation (C-5) 
are constant in our problem, 
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(a constant) 



(c-6) 



it 



FI 




FIGURE 0-1 

Heat Balance for an Element of Fluid 



Therefore, the temperature gradient in the direction of 
flow is a constant independent of radial location. 
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APPENDIX D 



NUSSELT NUMBER 



Nusselt number* Nu is a useful expression in comparing 
our problem results. It is function of the temperature grad- 
ient at the wall to a constant average temperature gradient^. 
Considering a heat balance on a unit length of tube; 



1 



i n 



rrp 

f 



_ q rrD ~ o u m c-pf\ rrb 

I, \a./ \ * I I 



L t 



I'W \ m 

where q w " is taken positive if heat is being removed^. Sub- 
stituting in equation (D-l) the equations 

Qv;" - 11 Oin 



(D~l) 



and 



gives 






m - 



/ 6 k 

Y'b 1 



11 _ hi) <L 31 

On j( f faUmtyR 



I 



(I)~2) 



Substituting the definition of F in equation (0-2) and rear- 
ranging terns gives 

LL 

(0-3) 

From this equation, as (J)^ approaches zero, Nu will ap- 
proach infinity; and when F -- 1 , Nusselt number will equal zero. 



Nu = hD/k = -4— (1-F) 
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APPENDIX E 



REYNOLDS NUMBER 



Reynold number Re is the ratio of the inertia forces to 
frictional forces and is a measure of turbulence . 

Considering a dynamic force balance on an element of 
fluid, the following expression can be developed: 



& - ~ x # fy- 



ix dx 



0 



tVn c j 



where -f- is for upward flow and - is for downward flow 
Prom equation (13) the pressure-drop parameter C is 



8 



c = 



(jf z fw I ) 



3 l y>\ u m 

and substituting equation (E~l) in equation (B-2) gives 



c = 






r*r 23? 
1_D_ 



V/ 



~^( " fm ) ~fc'1 



J 

It J 



m 

!/ 



3 2 A-i ^ 

Using equations 

Tw = ( «T 

fw — fryi — ^ 0 

On) - ^ D 2 $ m/ f <c ;( / 

and the definitions of Re, Re*, Ra, and P, we obtain 

c - *±S£L &Jm_ . 

8 /? e 3 £ F 

Solving for Re gives 

ho - LISsTt 

SC + Ra 

vzhich is equation (36) in the Analysis, 
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(E-2) 



(E-3) 



(E-4) 



(E-5) 



APPENDIX P 



VOLUME HEAT SOURCES 

The term q n * accounts for the volume heat sources in 
our problems. Since only uniform volume heat sources are 
treated, q MI is independent of direction. 

The nondimensiona'l volume heat source parameter is 



F - Pm C P Ur/) R (F-l) 

~ <^"i 

as defined by Hallmark. The physical significance of P can 
be described as the ratio of the thermal energy convected 
dovmstream, per unit volume to the heat generated in the fluid 
per unit volume. Considering a heat balance on a unit length 
of pipe v/ith constant v/all heat addition 



1 



I/I 



jr Q! 

L t 



c j" rrD ~0u m Cyj/}rr_D 

KW Xm ' Li 



(F-2) 



L t 



and introducing equation (P~l) gives 



P 




(P-3) 



The meaning of the various values of F can be seen from 
equation (P-3) . When q t " i: 0, v/hich corresponds to an insulated 
tube v/all, then F = 1. The walls should be hotter than the 
bulk of the fluid because heat is convected away from this 
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region less rapidly than in the center, and because of uniform 
volume heat source generation. Fluid velocities should there- 
fore tend to be higher in the vicinity of the wall than near 
the center line. 

When F < 1, it means that heat is being removed at the 
walls because more heat is being generated internally than is 
convected downstream. The velocity at the wall should be slow- 
er than the velocity at the center line in this case since heat 
is being removed and the density of the fluid is higher at the 
wall. 

When F>1, it means that heat is being added at the walls 
because heat is being convected downstream faster than it is 
being generated. The velocity in this case should rise faster 
near the walls than near the center line since it is hotter at 
the walls . 

F = 0 will occur whenever there is no net through-flow 
(u m r 0), no flow (u = 0) , or constant wall temperature (AzO) . 
An investigation is being conducted on the general problem 
with a constant wall temperature by Page 1 ?. 

P : oo would correspond to no volume heat sources. An 
investigation is being conducted on this general problem by 
Griffin 16 . 
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